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ABSTRACT 



We consider stochastic particle acceleration in plasmas around stellar mass 
black holes to explain the emissions above 1 MeV from Galactic black hole 
candidates. We show that for certain parameter regimes, electrons can overcome 
Coulomb losses and be accelerated beyond the thermal distribution to form 
a new population, whose distribution is broad and usually not a power law; 
the peak energy of the distribution is determined by the balance between 
acceleration and cooling, with particles piling up around it. Radiation by 
inverse Compton scattering off the thermal (from background) and non-thermal 
(produced by acceleration) particles can in principle explain the hard X-ray 
to gamma-ray emissions from black hole candidates. We present model fits of 
Cyg X-1 and GRO J0422 in 50 keV - 5 MeV region observed with OSSE and 
COMPTEL. 



Subject headings: acceleration of particles — black hole physics — radiation 
mechanisms: non-thermal — thermal — gamma rays: observations — theory 
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1. INTRODUCTION 

Most of Galactic black hole candidates (GBHCs) show X-ray spectra well fit by a 
thermal Comptonization model with temperatures ~ 50 - 100 keV and Thomson depths 
of a few (e.g., [Harmon et al. 1994| ; pijiang 1993|) . OSSE and COMPTEL experiments 



on Compton Gamma-Ray Observatory (CGRO), however, have recently revealed that 
persistent gamma-rays > MeVs are being produced in some GBHCs, in particular Cyg X-1 
( [Johnson et al. 1993| ; [McConnell et al. 1994a[ ) and GRO J0422+32 ( [van Dijk et al. 1995[ ). 



These gamma-ray tails are hard to fit with a single component Sunyaev & Titarchuk (1980) 
thermal model (or a more recent model by Titarchuk 1994). Hence non-thermal processes 
are strongly hinted by these observations. 

Models of non-thermal pairs have been studied with the emphasis on the power-law 
X-ray emission from active galactic nuclei (AGNs) (e.g., [Fabian et al. 1986[ ; [Light man 



Zdziarski 1987[ ; [Svensson 1987t jCoppi 1992| ), and they can be also applied to plasmas around 



stellar mass black holes. Those models assume that either mono-energetic or power-law 
leptons with a large Lorentz factor (7 ^ 1) are injected and they initiate cascade processes 
such as production and Compton scattering (see pvensson 1994[ for a recent review). 
However, those models did not specify any acceleration mechanism, thereby lacking the 
self-consistency in determining the particle distributions and photon spectra. 

Motivated by the MeV photons from GBHCs, we study the role of stochastic particle 
acceleration in accreting plasmas near GBHCs. The mechanisms of wave-particle resonant 
interactions that lead to particle acceleration have been directly observed in solar-wind 
(e.g., [Marsch 1991| ) and extensively studied in the context of solar flares (e.g., [Melrose 1974 ; 



Ramaty 1979[ ; [Hamilton fc Petrosian 199^ ; [Miller fc Roberts 1995[ ). Stochastic acceleration 



has been applied to diffusive shock acceleration and cosmic rays ( [Schlickeiser 199^ ), the 
lobes of radio galaxies (e.g., [Lacombe 1977^ jAchterberg 1979| ; [Eilek fc Henrikson 1984[ ), and 
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recently, the central regions of AGNs ( Dermer, Miller, fc Li 1995 ). Here in this Letter, we 
couple the particle energization and radiative processes to determine self-consistently the 
steady state particle and photon distributions in GBHC environments. We then present fits 
to gamma-ray data of GBHCs observed with OSSE and COMPTEL. 



2. THE MODEL 

As a first approximation, we average all plasma and magnetic field properties over an 
ad hoc spherical volume around the black hole of radius R. In this homogeneous spherical 
emission region, we treat the relevant radiation and particle acceleration/heating processes 
comprehensively, but ignore details of the accretion flow, primary energy generation and 
viscosity, etc., by assuming naively that the total GBHC luminosity is released uniformly 
throughout the spherical volume. This approach allows us to first concentrate on the 
microphysics. We assume that the plasma cloud is in a steady state, and solve the 
coupled steady-state kinetic equations of particles and photons (e.g., [Lightman fc Zdziarski 



1987|) . Our approach differs from previous studies in that we include the diffusion and 
systematic acceleration terms in the particle kinetic equations, both of which are derived 
from the stochastic acceleration. We also assume that the particles consist of a thermal and 
non-thermal component, coupled by Coulomb interactions. 

The evolution of the non-thermal particle distribution in momentum space, n^t{p), can 
be described by the Fokker-Planck equation as 

dn^t d 



dt dp 



d_ 

dp 



2 

-D{p)n^t + print 
p 



where p = |p|/meC = P'j, 7 is particle's Lorentz factor, and /? = (1 — 1/7^)^/^. Here, 
D{p) is the particle diffusion coefficient in momentum space, tgsc is the energy-dependent 
escape timescale, p is the cooling losses (both inverse Compton scattering and Coulomb 
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losses are included), Qinj(l') is the injection including sources from both the thermal 
background plasma and pairs production, and A{p) represents the removal of leptons 
due to pair annihilation (e"^ + e~ — > 7 + 7). The non-thermal particle density is given 
by / nnt{p)dp. A direct consequence of including D{p) (also the difference from previous 
study) is the development of a population of non-thermal leptons which is accelerated out 
of the background thermal plasma. However, below a certain momentum pthr, which we 
assume to correspond to the plasma temperature Te and 7thr = {Pthr + 1)^^^ = 1 + 46 (see 
e.g., Zdziarski et al. 1990 ), the particle distribution should be Maxwellian, possibly due to 
either waves are strongly damped at those energies or "thermalization" timescale (such as 
Coulomb collisions) becomes the shortest. Thus all the particles with p < pthr are grouped 
into the thermal population. We concentrate on the steady state solutions {d/dt = 0). 

We admit that the actual distribution of the turbulence are subject to great uncertainty. 
We assume that turbulent plasma waves are generated uniformly in the cloud and that there 
is a power-law distribution of turbulence energy density in wave number W{k) oc /c"'', where 
k = \k\, with k extending from kmin (which is taken to be ~ 1/i?) to the value near electron 
cyclotron resonance, so that it covers both the Alfven and whistler branches. Both D{p) 
and tesc can be evaluated from the conditions of accreting plasma (see recent treatment in 
Dermer, Miller, & Li, 1995). Here we combine the results from both the electron/whistler 
and electron/ Alfven interactions, and obtain the pitch-angle averaged momentum diffusion 
coefficient as 



TT 



1 PJ 



D(p) = -(g-l) — ^Cwave hr 



^dyn /5 



q-2 



RJ 



Ia{ko, ki) p'^ , 
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Ia{ko, ki] 
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with tdyn = -R/c, Cwave = Uw/Ub (the ratio of wave to total magnetic energy density), 
fL = meC^/eB (the electron gyroradius), ki = (mp/meY^'^ / Pa, ko = {mp/me)/p, and 



- 6 - 



c/^A = B / (ATmpmpY'''^ being the Alfven velocity, where Up is the proton number density. 
Here B is obtained by assuming equipartition between thermal gas pressure and magnetic 
pressure, and in most cases we study here. Pa ~ 0.01. 

The importance of stochastic acceleration can be illustrated in Figure 1, where we plot 
the various timescales (in units of tdyn = R/c) as functions of the dimensionless particle 
momentum p for a test electron passing through a plasma cloud with size i? 1.5 x 10^ 
cm, temperature G = kTf,/ {111^.0^) ^ 0.2, and varying density Hp. Here we assume that 
magnetic energy density Ub is in equipartition with thermal gas pressure (upQ) and 
Cwave = 0.1. Also, Mph = f/ph/f^B, whcrc f/ph is the soft photon energy density. All 
timescales are evaluated from energy change rates. The acceleration rate is calculated from 
{d'j/dt)i^c = [p'^PD(p)] and equation (^. The Compton+synchrotron cooling rate is 
I {d'y/dt)c+s\ = + 'Uphj'TpQp'^, where Tp = a^^UpR and ctt is the Thomson cross section. 

The cooling rate of lepton-lepton Coulomb interactions is taken from Dermer & Liang 
(1989). The diffusive escape timescale is typically the longest among all timescales and 
is not shown here. It is evident that, for certain parameter regimes, acceleration is more 
efficient than the cooling processes so that it must be taken into account in determining the 
particle distributions of emitting plasma. 

The energy flow rates for various components in the system are outlined as follows: (A) 
The total gravitational energy released by accretion is characterized by the compactness £, 
which is given by ^ = AtkjiR^u/ {"inif.c'), where u is the energy input rate per unit volume; 
(B) A fraction of gravitational energy, £nt = Cnt^? is used to stochastically accelerate some 
leptons through wave-particle resonant interactions; (C) The rest of the input energy, 
£th = (1 — Cnt)^, is used to heat up thermal leptons; (D) Soft photons, with a blackbody 
spectrum of temperature T^, are injected in the cloud uniformly with the compactness 4 
(the source of soft photons can be optically thick accretion disks). 
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In determining the particle distribution and emission spectra, we have included most 
of the relevant microphysical processes, all of which are put into a computer code (|Li 1995^ 



see also [Kusunose fc Mineshige 1995| and [Lightman fc Zdziarski 1987|) . These include (A) 



Non-thermal processes: Compton cooling. Coulomb collisions and annihilation between 
non-thermal and thermal leptons, the flow of non-thermal pairs becoming thermal pairs 
due to cooling, and the escape; (B) Thermal Processes: Coulomb collisions between protons 
and thermal e^, thermal bremsstrahlung and Compton scattering by thermal e^, and 
annihilation emission by thermal e^; Thermal Compton scattering is calculated using the 
Fokker-Planck equation (Kompaneets equation) with the diffusion coefficient given by 
Prasad et al. (1988) for isotropic radiation field. When Te = a^neR < 2 and Te > 100 
keV, however, we use a diffusion coefficient for the collimated radiation field, based on 
the argument by Titarchuk (1994). (C) Radiation Processes: Additional processes include 
the injection of soft photons is, absorption due to 7 + 7 interactions, and the escape from 
the system; (D) Pair Balance: Thermal pairs are assumed in pair balance, i.e., the pair 
production rate is equal to the annihilation rate for p < j^thr, and escape is neglected. 

To summarize our model, for given i, proton density rip, and R, distributions of steady 
state particles and photons for a plasma cloud with the processes described above are 
obtained for given parameters such as ent, is, Ts, and q. The temperature of the plasma 
Te is mainly determined by the balance among external heating (by protons) ith. Coulomb 
heating by non-thermal leptons, and cooling by thermal Compton scattering. The Thomson 
scattering depth for thermal leptons Te = a^UeR (usually different from Tp = a^UpR due 
to pairs) is also obtained. The level of the turbulent plasma waves Cwave is unambiguously 
determined from the relation £nt = (47rcrT-R^/3c) / dpnat{p){d'y /dt)s.cc, where {d'y/dt)acc 
is proportional to Cwave and n^tip) is obtained from equation (|l|). The advantage of this 
approach is that we do not need to artificially assume the injection of relativistic leptons as 
was done previously, because the content of non-thermal particles is now self-consistently 
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determined from the stochastic particle acceleration. 



3. RESULTS 

Figure ^ summarizes our main results. Photon spectra of Cyg X-1 and GRO J0422 
from OSSE (cross, |Fhlips et al. 1995| , |Kroeger et al. 1995| ) and COMPTEL (filled square, 



McConnell et al. 1994b|, [van Dijk et al. 199 51) are shown in the left panels. 



(1) Cyg X-1. Our best-fit model spectrum (left) to Cyg X-1 and the corresponding 
particle distribution (right) are shown as the solid curves in the upper panel of Figure 
0. Other parameters are given in Table |I[ Clearly a deviation from Maxwellian of the 
particle distribution occurs at high energy tail. This occurs because the "upward" motion 
(in momentum space) by acceleration is balanced with the "downward" motion by cooling 
at a specific value of momentum, at which particles tend to "pile up" , forming a "bump" 
in the particle distribution (see also Schlickeiser 1984). It is from this "bump" component 
that most of the emissions above 1 MeV are produced via inverse Compton scattering. 

Compared with the thermal models, only i^t and q are the new additional parameters. 
While the dependence of the results on q is rather weak, we demonstrate the effects of 
varying i^t in the upper panels of Figure |^ (dotted and dashed curves). Dotted curve 
(^nt/^ = 10~^) indicates that photon spectrum from Cyg X-1 can not be fitted by pure 
thermal emissions. As i^t increases, more pairs are produced (dashed curve in particle 
distributions), so is > MeV emission in this case (but see below). Note that the pair 
fraction (1 — Tp/Te) in this case is small (~ 3.4%). 

(2) GRO J0422. A good fitting to GRO J0422 is shown by the solid curves in lower panel 
of Figure |^. For M = IOMq, we find that it only needs int/^ = 0.06, and the particle 
distribution is dominated by the thermal component with only a slight deviation from the 
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Maxwellian distribution at higher energy due to acceleration. This is also evident from 
the fact that MeV emission from J0422 is much weaker, compared to Cyg X-1. We have 
also tried to fit J0422 by choosing a different mass but keeping other parameters to be 
the same as for Cyg X-1 (as suggested by the referee). Dashed and dotted curves in the 
lower panel show other fittings (though not as good as the M = IOMq case) to J0422 with 
M = 6MQ,i= 15 and M = 4:Mq,£ = 22.5 (see Table |ip. Note that the source luminosity 
is (X i M for fixed R so it is the same for all trials. The particle distributions are now 
essentially thermal. Again, other parameters are given in Table 0. 

Figure ^ also indicates an important feature of our model, i.e. higher £nt/^ does not 
necessarily result in higher fiux of gamma-ray emissions. In the case of J0422, the high 
compactness (smaller mass) causes the cooling becoming increasingly important at high 
energies, destroying the otherwise suprathermal population and merging them into thermal 
bath. Notice that the cutoff energy in photon spectra decreases and the number of pairs 
increases dramatically (~ 26% for 4Mq) as mass decreases. 

4. CONCLUSIONS AND DISCUSSIONS 

We have examined stochastic particle acceleration via wave-particle resonant 
interactions near Galactic black holes. We find that under certain conditions, stochastic 
electron acceleration can overcome both Coulomb and Compton losses, resulting in a 
suprathermal population. Preliminary model spectra show good fits to the recent OSSE and 
COMPTEL observations of Cyg X-1 and GRO J0422. We find that Cyg X-1 has a much 
higher component in gamma-rays than J0422, as it is also evident in the data. This can be 
understood in terms of the higher £nt/^ and lower compactness of Cyg X-1, compared to 
J0422, if they have the same black hole mass. On the other hand, the high energy emissions 
from J0422 can be fitted by different masses of the putative black holes. But these fittings 
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are different from the standard thermal (i.e. e-p plasma only) models because a large 
amount of pairs are produced and thermalized due to the high compactness. Since both 
sources are highly variable, especially in gamma-ray emissions, we expect the non-thermal 
energy content £nt/^ might vary but the detailed microphysics of how the system partitions 
the energy flows is presently unclear. 

The non-thermal population we obtain here is relatively soft with few pairs, not 
capable of explaining either the transient MeV bump discovered by Ling et al. (1987) or 
the COMPTEL power-law tail if it extends to much higher energies. This softness is due 
to our assumptions of a copious soft photon source permeating the entire emission region 
and the coexistence of thermal and non-thermal particles in a single homogeneous volume. 
Coulomb interaction efficiently transfers the non-thermal energy to the thermal reservoir 
which strongly limits the existence and energy content of non-thermal particle population. 
To achieve much harder non-thermal populations and emissions we need to postulate the 
existence of physically distinct regions for the thermal and non-thermal components, either 
in the form of a thermal disk plus non-thermal corona or jets, or thermal outer disk and 
non-thermal inner torus. We plan to undertake such ventures in future work. But it 
is highly encouraging that a simple-minded single-component spherical model with few 
parameters is already capable of explaining from flrst principles the quiescent (low-hard 
state) emissions of Cyg X-1 and GRO J0422 detected by CGRO. 

We note, of course, that there are many other alternative interpretations of the hard 
tail, including the pion decay model (|Jourdain fc Rogues 199^ ) and pair-dominated hot 



cloud model invoking quenching of the soft photon source ( [Liang fc Dermer 1988| ). It is 



possible that particles accelerated by shocks may also account for the > MeV emissions. 
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Fig. 1. — Timescales of acceleration {solid), Coulomb loss (dotted), and Compton loss 
(dashed) as functions of particle's dimensionless momentum p = 7/?. The scattering depth 
and soft photon energy density, indicated respectively by Tp and tiph (see text for definitions) 
for a given size R, are varied as (rp, u^h) — (0.05, 1), (0.5, 1), (0.5, 10), and (2, 0.1), for plots 
(a), (b), (c) and (d), respectively. Sufficient acceleration occurs in plots (a) and (b), but 
Compton/syn. cooling and Coulomb cooling prevent acceleration in (c) and (d), respectively. 

Fig. 2. — Upper panels: Model fit of the emission spectrum of Cyg X-1 (left) and the 
corresponding particle distribution as function of dimensionless momentum p (right). Data 
are from OSSE and COMPTEL (squares) experiments. Sohd curve (int/^ = 0.15) fits the 
data. Dotted and dashed curves are for £nt/^ = 10^^ and 0.5, respectively. They depict the 
effects of changing ent only, with other parameters fixed. Particle distributions for higher 
^nt/^ (solid and dashed curves) show deviations from Maxwellian distributions in the high 
energy tail due to acceleration. Lower panels: Same as the upper panels but for GRO J0422. 
Again, data are from OSSE and COMPTEL. The solid, dashed and dotted curves are for 
M/Mq — 10, 6, and 4, respectively. Sohd curve fits the data best with int/'^ — 0.06, and 
the corresponding particle distribution is very close to Maxwellian. Particle distributions for 
M/Mq = 6, 4 are essentially thermal and a large amount of pairs are produced (dashed 
and dotted curves). 
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Table 1. Model Parameters and Results 



Object 


dist. 


mass 


I 


^/4dd 






Ts 


Tp 


q 


R 




Te 


Cwavc 


Cyg X-l'' 


2.5 


10 


4.5 


0.02 


0.15 


0.07 


0.1 


0.7 


5/3 


80 


139 


0.725 


0.059 


J0422" 


2.5 


10 


9.0 


0.03 


0.06 


0.08 


0.1 


0.7 


5/3 


80 


134 


0.741 


0.065 


J0422^ 


2.5 


6 


15 


0.05 


0.15 


0.07 


0.1 


0.7 


5/3 


80 


119 


0.869 


0.08 


J0422'= 


2.5 


4 


22.5 


0.08 


0.15 


0.07 


0.1 


0.7 


5/3 


80 


111 


0.944 


0.07 


''distance in kpc 


mass 


in Mq, temperature in 


keV, R in 


units of GM jf? 











''parameters for the solid curve for Cyg X-1 in Figure |^ 

'^''''''parameters for the sohd, dashed and dotted curves for J0422 in Figure ||, respectively 



